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Neutrino diffusion in the pasta phase matter within the Thomas-Fermi approach 
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The behavior and properties of neutrinos in non-uniform nuclear matter, surrounded by electrons 
and other neutrinos are studied. The nuclear matter itself is modeled by the non-linear Walecka 
model, where the so-called nuclear pasta phase is described using the Thomas-Fermi approximation, 
solved in a Wigner-Seitz cell. We obtain the total cross-section and mean-free path for the neutrinos, 
taking into account scattering and neutrino absorption, and compare the final results for two known 
kind of model parametrizations: one in which non-linear effects in the strong sector are explicitly 
written in the model Lagrangian and another one in which the coupling constants are density 
dependent. The solution for this problem is important for the understanding of neutrino diffusion 
in a newly born neutron star after a supernova explosion. 
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1. INTRODUCTION 

Neutrinos are elementary particles that interact with 
other particles only through the weak force, which makes 
its scattering by matter very unlikely. On the other hand, 
this pure weak force scattering can reveal some aspects 
of the structure of matter that other stronger interac¬ 
tions can not. Also, this feature represents an advan¬ 
tage in the sense that neutrinos carry information about 
the target along great distances, as is the case of neu¬ 
trinos produced in the interior of stars that reach earth 
detectors. For those reasons, the study of neutrino in¬ 
teraction with matter is at the same time very fruitful 
but very challenging. In what refers to the scattering 
of neutrinos by hadronic matter, despite of the experi¬ 
mental as well as theoretical problems, some important 
progress has been made as in the case of the Karmen col¬ 
laboration [i] and Los Alamos results [ 2 } . More recently, 
long-baseline experiments are under way, like MiniBoone, 
NOvA and other similar accelerator experiments [|]. Al¬ 
though the main purpose is to obtain further information 
on Standard Models and on neutrino oscillations, a pre¬ 
cise knowledge of the neutrino-hadron interaction as well 
as the hadronic structure of the targets 0] is necessary. 

Another important source of information comes from 
the physics of the supernova core collapse and the evolu¬ 
tion of a neutron star. The important ingredient in this 
case is the propagation of neutrinos in such a medium, 
which contains not only baryons but other leptons. In 
particular, for baryonic densities below the nuclear satu¬ 
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ration value, structures known as pasta phases [H| are ex¬ 
pected, which in some sense resemble the baryonic struc¬ 
ture in nuclei, but now embedded in a neutron gas and 
an electronic distribution. Actually, a delicate competi¬ 
tion between Coulomb and surface energies determine the 
most favorable final inhomogeneous structure. The pres¬ 
ence of such structures presumably have a non-negligible 
role in the neutrino diffusion which is a key ingredient 
for the modeling and simulation of the core collapse su¬ 
pernovae mechanisms [f§-[8[ ■ 

One way to obtain the pasta phase is to solve the 
problem considering charge neutral Wigner-Seitz cells 
of appropriate geometries containing neutrons, protons 
and electrons within a variational approach, using both 
relativistic and non-relativistic mean-field calculations. 
Most of the recent applications in this case have used 
the Thomas-Fermi approximation mm , but Hartree- 
Fock calculations fill and three dimensional Thomas- 
Fermi calculations jl2[ can also be found in the recent 
literature, all within the Wig ner-Seitz approximation. 
Recently, it was shown in [13| that performing a pasta 
calculation taking a large enough cell to include several 
units of the pasta structures, different distributions of 
matter from the usual ones considered within the Wigner 
Seitz (WS) approximation could be energetically favored 
in some density ranges. Another approach to this prob¬ 
lem that also goes beyond WS approximation is based in 
the so-called quantum molecular dynamics EMI. 

Here we follow the Thomas-Fermi results as described 
in 0 in order to generate the pasta-phase structure, 
where a relativistic model lagrangian is the starting 
point. The self-consistent calculation is performed con¬ 
sidering matter in /3-equilibrium, where just protons, 
neutrons, electrons and neutrinos are present. The total 
neutrino cross-section for each kind of particle is calcu¬ 
lated, as a function of the density, taken the correspond¬ 
ing geometry for the considered density, i.e., droplet, 
rod, slab, tube or bubble. The neutrino mean-free path 
(NMFP) energy and temperature dependence for selected 
values of the density is also presented. 

The neutrino cross-section calculation includes the ef¬ 
fects of strong interaction, accounted for in-medium mass 
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and energy shifts and degeneracy effects, based on the 
formalism previously developed for homogeneous nuclear 
matter [ 28 j, with the difference that now the strong 
and electromagnetic potentials as well as the nucleon 
masses and energies are position dependent. Conse¬ 
quently, the result is similar to the free space transi¬ 
tion amplitude, and the uniform matter result becomes, 
straightforwardly, a particular case of our expressions for 
the cross-section. 

We will not discuss in the present paper the impor¬ 
tant sources of suppression and enhancement due to in¬ 
medium correlation 0. Our main objective is to un¬ 
derstand how inhomogeneous matter affects the neutrino 
cross section through the in-medium effects. 

As in [b§], we have considered elastic scattering and 
neutrino absorption in our derivation. In other words, 
neutral current and charge changing processes are in¬ 
cluded and their relative importance to the total cross- 
section is discussed. In what follows, an outline of the 
formalism is presented in section[2] numerical results and 
discussion are provided in section 0 and the conclusions 
are in the final section 2) The details of the pasta phase 
calculation can be found in the cited references and some 
details of the cross-section calculation are shown in the 
Appendix. 


2. FORMALISM AND MODEL 
PARAMETRIZATION 

We start with a model Lagrangian density that in¬ 
cludes electrons, neutrinos, nucleons, the sigma, omega, 
rho and delta meson fields and the electromagnetic inter¬ 
action, given by EM: 


C— 'y ' Ct +C (T +C ul -\-Cp+C U] p-\-Cs-\-C~ 1 -\-C e ( 1 ) 

i=p,n 

where the nucleon Lagrangian reads 

( 2 ) 

with 

iD " = id^ — Y v V^ — — e—^-^-A^, (3) 

m* = M-r s <f>-r s r-8. (4) 

The meson and electromagnetic Lagrangian densities 


are 

= \ - m 2 s cj ) 2 - ^ 

A, = ^-isi F r + m5ij 1 F + ^(i;F)2j 

Cp = + 

Cs = \{d p 8d»8-m 2 5 8 2 ) 

r — _i Jr 

L up = A^X-b^C^)), 

where = d p V v -d v V u , , B pv = dph v -d v hp,-T p {x 
by) and F pv = d p A v — d v A p . The electromagnetic cou¬ 
pling constant is given by e = y47r/137 and r is the 
isospin operator. Finally, the electron and neutrino La¬ 
grangian densities read: 

£ e = V>e [lp(id P + e A P ) - m e \ i>e 
Cv = tV 

The weak boson contributions to the above La- 
grangians, affect the solutions of the corresponding Euler- 
Lagrange equations, in a completely negligible way, due 
to their huge masses and relatively small energy range, 
and are only relevant in order to obtain the neutrino cross 
sections, so they are not included in the above equations 
in the determination of the equations of state. 

The solution of the corresponding Euler-Lagrange 
equations are obtained in the mean field self-consistent 
Thomas-Fermi (TF) approximation, as explained in ref¬ 
erences E3, 0. We have considered two types of 
parametrizations: one in which the couplings are taken 
to be constants (NL) and another one in which they are 
density dependent (DD). In the last case, the terms pro¬ 
portional to k, A, £ and A are set equal to zero. Within 
TF the system is considered locally uniform and the main 
output are the densities, which are position dependent. 
Explicitly, we have for the baryonic, scalar, isoscalar, 
scalar-isoscalar, electron and neutrino densities: 

P( r ) = P P ( r ) + Pn(r) = ; 


Ps (r) = p Sv (r) + p Sn (r) = ; 


Pair) = p P (r) - p n (r) = ; 


Ps3(r) = p S p{r) - p sn {r) = ; 
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Pe(r) = (i>ti> e y, 

Pu{r) = ; 

with 

Pi(r) = (2^)3 / d3fc (%*( T ) “ Vki( T )), i=p,n,e,v; 

( 5 ) 

^( r ) = (^3 J J^riVkiiT) + fjki(T)), i=p,n ; 

( 6 ) 

where E* = \/k 2 + M* 2 , k is the momentum and 7 = 2 
is the spin multiplicity. For a given temperature T, the 
distributions are: 

Vki(T ) = {exp [(£* - /+)/T] + l} -1 ; 

Vki{T) = {exp [(A - /xO/T] + l}" 1 ; (7) 

with Ei (Ei) and /+ (pi) being the particle (antiparticle) 
energy and chemical potential, respectively. The particle 
(antiparticle) energy depends on the mesonic fields and 
is position dependent. 

Once the densities are determined we proceed for the 
calculation of the total neutrino cross-section (a). We 
follow here the procedure discussed in [l 8 | . For a collision 
1 + 2 —> 3 + 4 we may write: 


f d 3 k 2 

CO 

CO 

'TS 

f d 3 ki (27t) 4 

) (2t r) 3 j 

f (2 n) 3 J 

(27t) 3 \v-i - v 2 \ 


S\Pi +P 2 -P 3 - PMT)(1 - jfc(T))( 1 - 774(T)). 

{(V + A) 2 (l - v 2 cos(0 12 ))(1 - U 4 COs( 6 » 34 )) 

+ (V — A) 2 (l — V4,cos(0\i))(l - V2 COs(92 3 )) 

M* M* 

-^±(V 2 -A 2 )(l-cos(9 13 ))}. ( 8 ) 

E 2^4 

In the above expression Pi is the particle four mo¬ 
mentum, Vi = \ki\/E*, Gf is the Fermi constant 
and V, A are the weak current vector and axial vec¬ 
tor couplings, respectively, and depend on the target 
particle and on the exchanged weak boson. For the 
neutrino-electron (neutrino) cross section we replace M* 
by TO e (zero). The factors (1 — 773 (T)) and (1 — rn(T)) 
are due to final state Pauli blocking effects. We have 
considered here scattering through neutral current and 
charge-changing processes as well as neutrino absorption 
by the neutrons. The explicit expressions and definitions 
for each case are shown in the Appendix. Through the 
analysis of those expressions we may conclude that the 
integrand in the cross section is r dependent, once the 


potentials and the effective mass are position dependent. 
Note that the same expression can be used for the calcu¬ 
lation of the cross section in the infinite nuclear matter, 
for which the potentials and effective mass are not r de¬ 
pendent. 

The nucleon internal structure was taken into account 
in our calculation, just multiplying the couplings V and 
A by the vector and axial-vector nucleon form factors as 
explained, for instance, in | 2 l| . 

The neutrino mean-free-path is then obtained consid¬ 
ering the cross-section for a Wigner-Seitz cell divided by 
its volume 



3. RESULTS AND DISCUSSION 

In what follows we show our main results for the pasta 
phase, using both Lagrangian parametrizations with con¬ 
stant and density dependent couplings, as explained be¬ 
fore. For the parametrizations with non-linear terms in 
the mesonic sector and constant couplings, we choose the 
GM3 [U, NL3 mi, NL3wp H and the FSUGold jUJ 
sets of parameters. The sets NL3 and NL3u;p only differ 
on the density dependence of the symmetry energy and 
will allow the discussion of the effect of this quantity on 
the neutrino mean-free path. For the density dependent 
case we take the TW [Jfj] and the van Dalen et al. [27j 
sets. The last one was recently taken for neutrino mean- 
free-path studies Et Table [I] shows the properties of 
the models for symmetric nuclear matter at saturation 
density p 3 and zero temperature, and in Fig. [T|we have 
plotted the symmetry energy of the models versus the 
density, a quantity that will be necessary to discuss the 
results in the following. Models Dalen, GM3 and NL3 
have a smaller symmetry energy in all or part of the den¬ 
sity range below 0.1 fm -3 . A smaller symmetry energy 
allows the system to have a larger isospin asymmetry. As 
we will see next these three models are precisely the ones 
with smaller proton fractions in the pasta phase matter, 
and this will induce a noticeable effect on the neutrino 
mean-free-path. 

The free energy per particle obtained with the 
FSUGold parametrization for /^-equilibrium matter with 
trapped neutrinos for a fixed fraction of leptons Yj, = 0.4 
is shown in Fig. [2j where the different geometries are 
identified by different colors. The uniform matter result 
is also shown, and, as expected, has a larger free energy 
than the pasta phases. Although the transition densities 
between geometries may differ slightly, depending on the 
parametrization used [28|, the main behavior does not 
change significantly compared to the case shown. 

Figures 13a) and H] display, respectively, the proton 
fraction and the chemical potential for all particles in¬ 
volved using again the FSUGold parametrization to de¬ 
scribe neutrino trapped /^-equilibrium matter. All these 
results were obtained for a temperature T = IMeV. One 
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TABLE 1: Properties of the models parametrizations for in¬ 
finite symmetric nuclear matter at zero temperature at satu¬ 
ration density: the saturation density po, the binding energy 
Eb/A, the incompressibility K, the nucleon effective mass 
M *, the symmetry energy a sy m and its slope L. 



po Eb/A K 

(fm^ 3 ) (MeV) (MeV) 

M*/M 

O'sym L 

(MeV) (MeV) 

FSUGold 

0.148 

16.299 

271.76 

0.6 

37.4 

60.4 

GM3 

0.153 

16.3 

240 

0.78 

32.5 

89.66 

NL3 

0.148 

16.3 

272 

0.6 

37.4 

118.3 

NL3tap 







CO 

o 

II 

s> 

0.148 

16.3 

272 

0.6 

31.7 

55.2 

Dalen 

0.178 

16.25 

337 

0.68 

32.11 

57 

TW 

0.153 

16.247 

240 

0.555 

33.39 

55.3 



FIG. 1: Symmetry energy as a function of the total baryonic 
density for all the models used in the present work. 

important effect of the clusterization of matter is the in¬ 
crease of the mean proton fraction. An effect of the larger 
proton fraction for matter with a fixed lepton fraction is 
the reduction of the neutrino fraction and, consequently, 
the neutrino chemical potential as seen in Fig. H) At low 
temperatures, when degeneracy effects are important, a 
smaller neutrino chemical potential will give rise to larger 
mean-free paths, and neutrinos will diffuse out of the star 
more easily. In Fig. [3Jb) the pasta phase matter pro¬ 
ton fraction is plotted for all the models studied. These 
results reflect the density dependence of the symmetry 
energy of the respective models: a smaller proton frac¬ 
tion occurs for the models with smaller symmetry energy 
at the baryonic density range of interest. In particular 
the three models with smaller proton fraction are GM3, 
Dalen and NL3, the three models that have the smaller 
symmetry energy for 0.05 < p < 0.1 fm -3 . 

The individual mean-free-path contribution for each 
particle type is shown in Fig. 0 as a function of ps, 
for the FSUGold and T = IMev. The curves labeled 



FIG. 2: Free energy as a function of the total baryonic den¬ 
sity using the FSUGold parametrization for uniform matter 
(black) and the pasta phase matter within droplet (red), rod 
(green), slab (blue), tube (purple) and bubble (light blue) ge¬ 
ometries. The temperature was taken as T = 1 MeV. M is 
the free nucleon mass. 

proton, neutron and neutrino are the contributions for 
elastic neutral current scattering, while the curve labeled 
electron has also an elastic contribution from charged 
current scattering. It is clear the dominance of the ab¬ 
sorption process for the total cross section. 

The sensitivity of the total mean-free-path to the 
parametrization used is presented in Fig. [G]together with 
the uniform matter result obtained with the FSUGold 
parameter set. At low densities the models do not dif¬ 
fer much. This is expected since all models have similar 
behaviors, in particular, the pro ton fraction is almost 
the same in all of them, see [231. The differences oc¬ 
cur precisely after the onset of the non-spherical pasta 
structures, ps > 0.02 fm - , which is a range of densities 
that is sensitive to the density dependence of the sym¬ 
metry energy. GM3 and Dalen have the lowest mean-free 
paths: a smaller proton fraction (see Fig. [3]) corresponds 
to a smaller electron fraction, which increases the neu¬ 
trino fraction and the neutrino chemical potential. All 
these facts favor the absortion process, and consequently 
the NMFP is smaller. Comparing NL3 and NL3u;p we 
also conclude that the softer density dependence of the 
symmetry energy gives rise to larger mean-free paths. 

Considering the scattering of the neutrinos from the 
individual nucleons of the pasta, as done in the present 
calculation, we see that the mean-free path increases a lot 
comparatively to the homogeneous matter case, mainly 
at low densities, when the clusterized matter presents 
much larger proton fractions, and, therefore, the ab¬ 
sortion process is less likely to occur. Other reason is that 
the processes are strongly suppressed by Pauli blocking 
effects inside the cluster, where the Fermi energy is larger, 
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FIG. 3: Proton fraction as a function of the total baryonic 
density, (a) using the FSUGold parametrization for uniform 
matter (black) and the pasta phase matter within droplet 
(red), rod (green), slab (blue), tube (purple) and bubble (light 
blue) geometries, (b) for all the models and the pasta phase 
matter. The temperature was taken as T = 1 MeV. 


while outside the cluster the space is almost empty. As 
expected, the larger the neutron and proton densities of 
the background gas the closer get the pasta and homo¬ 
geneous matter mean-free paths. 

The dependence of the results with the temperature 
can be seen in Fig. [3 Temperature has a strong ef¬ 
fect on the pasta structures which start to melt. Al¬ 
though within a TF calculation pasta structures still ex¬ 
ist at T > 10 MeV , according to [29|, if thermal fluc¬ 
tuations are considered the Wigner-Seitz cell structure 
is supposed to melt for T > 7 MeV. Temperature in¬ 
creases drastically the background gas of dripped par¬ 
ticles in the Wigner-Seitz cells and, therefore, the larger 
the temperature the closer come the pasta mean-free path 



FIG. 4: Chemical potential as a function of the total bary¬ 
onic density using the FSUGold parametrization for the pasta 
phase (points) compared to homogeneous matter (lines) for 
the neutrons (green), protons (red), electrons (blue) and neu¬ 
trinos (black). The temperature was taken as T = 1 MeV. 



FIG. 5: Mean-free-path as a function of the total baryonic 
density using the FSUGold parametrization for the pasta 
phase matter. Individual contributions for the electrons (red), 
neutrinos (green), protons (blue) and neutrons (purple). The 
absorption contribution is shown in light blue and the total 
value is shown in black. The temperature was taken as T = 1 
MeV and E v = /x„. 

to the homogeneous matter ones. The reduction of the 
mean-free path with the increase of the temperature is 
also expected, because the Fermi-Dirac distributions are 
smoothed when the temperature increases, making pos¬ 
sible more reactions. 

In Fig. [8] the NMFP is plotted as a function of 
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FIG. 6: Mean-free-path as a function of the total bary- 
onic density comparing various parametrizations for the pasta 
phase matter. TW (red), Dalen (green), GM3 (blue), NL3 
(pink), NLojp (grey) and FSUGold (orange). The FSUGold 
for homogeneous matter (black), is also shown. The temper¬ 
ature was taken as T = 1 MeV and E u = 




E v (MeV) 


FIG. 8: Neutrino mean-free-path as a function of the neutrino 
incident energy for the pasta phase with ps — 0.018fm _1 , 
comparing the parametrizations TW(red), Dalen(green), 
GM3(blue), NL3(pink), NLojp(grey), FSUGold(orange) and 
FSUGold homogeneous matter(black). The temperature was 
taken as T = 1 MeV (upper panel) and T = 5 MeV (lower 
panel). The triangles indicate the NMFP at the neutrino 
chemical potential, p„, of the respective model. 


FIG. 7: The neutrino mean-free-path as a function of the to¬ 
tal baryonic density within FSUGold parametrization for the 
pasta phase (points) compared to homogeneous matter (lines) 
for T = 1 MeV(red), T = 3 MeV(green), T = 5 MeV(blue) 
and T = 8 MeV(black) and E v = p v . 


the neutrino incident energy for the pasta phase with 
Pb = 0.018fm _ . We also include the homogeneous mat¬ 
ter result obtained within FSUGold. With triangles, we 
select the values of A at the neutrino chemical potential, 
at the two temperatures shown, T = 1 and 5 MeV. We 
expect that the neutrinos with energy around the chem¬ 


ical potential will provide the main contribution to the 
dispersion of energy in the system, because at low tem¬ 
peratures the system is practically degenerate and Pauli 
blocking factors will weaken other contributions. The 
highly degenerate regime is expected for pi/T 1, and 
this condition is essentially true for temperatures below 
the pasta melting temperature. The differences between 
the two temperatures is mainly explained by the drip of 
nucleons out of the clusters increasing the background 
gas, and reducing the differences between the pasta and 
homogeneous matter results. Temperature also decreases 
the mean-free path by almost one order of magnitude 
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E v (MeV) 


FIG. 9: Neutrino mean-free-path as a function of the neu¬ 
trino incident energy for the pasta phase with pB = O.OSfrrU 1 
and T = 1 MeV, comparing the parametrizations TW(red), 
Dalen(green), GM3(blue), NL3(pink), NLcjp(grey), FSUG- 
old(orange) and FSUGold homogeneous matter (black). The 
triangles indicate the NMFP at the neutrino chemical poten¬ 
tial, p v , of the respective model. 


both for the pasta and the homogeneous matter calcula¬ 
tion and this is explained by the opening of new transi¬ 
tions that are Pauli blocked at T close to zero. 

In Fig. [9]we consider a larger density corresponding to 
the slab geometry at T = 1 MeV. For these densities the 
gas of dripped neutrons is denser and the proton fraction 
smaller (see the blue region in Fig. [3] in comparison with 
the red region). Consequently, the NMFP decreases and 
comes closer to the homogeneous matter result. 


4. CONCLUSIONS 

We have studied the effect of the pasta phase, occur¬ 
ring in the inner crust of a neutron star, on the neu¬ 
trino mean-free-path (NMFP). The pasta phases have 
been obtained within a self-consistent Thomas-Fermi ap¬ 
proximation, as explained in references HI HI- Several 
relativistic nuclear models have been used to describe the 
pasta phase, both with nonlinear meson terms and con¬ 
stant couplings , and with density-dependent coupling 
constants. In particular, we were interested in discussing 
whether the properties of the models, such as the density 
dependence of the symmetry energy, would have some in¬ 
fluence on the NMFP. We have also studied the effect of 
the temperature. It should be stressed, however, that the 
present work is restricted to the density and the temper¬ 
ature range for which the pasta phases exist, e.g. T < 8 
MeV and ~ 0.0002fm -3 < p < p t where the density at 
the crust-core transition p t ~ po/2. We have considered 


both charged and neutral current reactions. It has been 
shown that the absortion process dominates the NMFP, 
but other processes can not be discarded. 

At low density, ps < 0.02fm -3 , where the pasta phases 
obtained within the different parametrizations have simi¬ 
lar properties [H}, the NMFP has a small dependence on 
the parametrization. As the density increases, namely at 
the layers close to the upper boarder of the inner crust, 
the dependence on the parametrization becomes more 
important. With a larger proton fraction, the absortion 
process is less likely to occur. Consequently, models with 
larger proton fraction result in larger NMFP, and since 
the proton fraction is closely related to the density depen¬ 
dence of the symmetry energy, its effect on the NMFP is 
not negligible. We have shown that the NMFP is larger 
in the presence of pasta phases than considering homo¬ 
geneous /3-equilibrium matter at the same density. Be¬ 
sides the fact that the pasta phases have a larger proton 
fraction, the processes are strongly suppressed by Pauli 
blocking effects inside the cluster, where the Fermi en¬ 
ergy is larger, while outside the cluster the space is al¬ 
most empty. When the density increases, the closer get 
the pasta and homogeneous matter mean-free paths. 

As the temperature increases, the pasta phase starts 
to melt and the pasta phase EOS approaches the infinity 
nuclear matter EOS, and, therefore, the NMFP in both 
cases get closer. At the temperatures we have considered 
the NMFP neutrinos are degenerate, and, therefore, the 
Pauli blocking factors ensure that only neutrinos close to 
the Fermi energy are involved at the dispersion of energy 
in the system. Comparing with the homogeneous case, 
we expect that neutrinos with less energy are involved 
in the presence of pasta phases, which give rise to larger 
NMFP. 

The NMFP as function of the neutrino energy for the 
different models do not differ much at low densities and 
temperatures, because the models have similar properties 
in this region of phase space. With the increase of the 
temperature and density, the dependence on the mod¬ 
els starts to appear. The great variation of the NMFP 
with the neutrino energy around the neutrino chemical 
potential is attributed to Pauli blocking effects. 

Our results imply that the effects of the pasta phase 
can not be neglected when calculating the NMFP at 
baryonic densities below O.lfm - and temperatures be¬ 
low 10 MeV. The parametrization used to described the 
pasta phase has also influence on the results, in particular 
due to the density dependence of the symmetry energy. 

We have also not considered coherent scatteringof the 
neutrinos from the pasta clusters as done in [6j, |7l JJJ . 
These processes are important when all neutrons respond 
coherently, for scattering by neutrinos with a wavelength 
of the size of the cluster, or energies E v < 75 MeV. The 
scattering processes discussed in the present work become 
important when the inner constitution of the clusters is 
distinguished by the neutrinos, and, therefore, for neu¬ 
trino with an energy above that value. 
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5. APPENDIX 


D = (ki + k 2 ) 2 + Ml- 
C — k\ T El — g p bo ; 

F = k\cos(9 3 ) + k 2 cos(9 2 3)- (12) 

The constants are: 

neutrino-proton: cy = 1/2 — 2 sen 2 (0 w ); ca = 1.23/2 
neutrino-neutron: cy = —1/2; ca = —1,23/2 
neutrino-electron: cy = 1/2 + 2 sen 2 (9 w ); ca = 1/2 
neutrino-neutrino: cy = \[2\ ca = V% 
absorption : gy = C; gA = —1.23 C 


According to equation (jS]|, with k\ = k\z and after the 
fc 4 integration in the scattering case: 
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|ui - v 2 \ 

k\ + El — k 3 

k\(l — cos{0 3 )) + El — k 2 cos{9 2 3) 
{(c v + C A ) 2 ( 1 - V 2 COs(0 2 ))( 1 - / 123 ) 

+ (Cy - C A ) 2 { 1 - V 2 COs(9 23 )(l - 3123 ) 
Ml Ml . 2 2 . . 1 
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ki + k 2 cos(9 2 ) — k 3 cos(9 3 ) 
E* ' 


_ fci ( El - k 2 cos(92)) 

3 k\ + El — k\cos(9 3 ) — k2COs(9 23 ) ’ 

For the absorption, we find: 

4 p p p tv p2n 

Cabs = J ^5 J d3r J J sen{9 3 )d9 3 J d(j> 3 kl- 

7? 2 (T)(1-?73(T))(1-7?4(T)) 

|«i - W 2 I 

^ tti 2 + A : 2 (C \/m 2 + fcf) 

- F^ml + k 2 

{(gv + 3 a ) 2 (1 - v 2 cos(9 2 ))( 1 - / 123 ) 

+ (9V ~ 9a) 2 { 1 - ffl23)(l - V2C0S(9 23 )) 

-^±(g 2 v -g 2 A )(l-COs(6 3 )} ( 11 ) 

where now 

__ 2 F(C 2 +m 2 e -D) 

3 4(kicos(9 3 ) + k 2 cos(9 2 3)) 2 — 4C 2 

yJ16m 2 e C 2 F 2 + 4 C 2 (C 2 + m 2 - D) 2 - 16m 2 C 4 _ 
4F 2 - 4C 2 ’ 


with C = 0.973 and sen 2 {9 w ) = 0.230. The collisions 
neutrino-electron and neutrino-neutrino can be repre¬ 
sented by two (first-order) Feynman diagrams each one, 
in such a way that the values of cy and ca can be re¬ 
defined to accommodate the different cross section con¬ 
tributions in a single expression, as given above. Also, 
we define: 


E* = Vk 2 +M* 2 . 

For the particle energy we have for the nucleon: 

Ei = E* + r„P 0 (r) ± \r p b 0 (r) + eA 0 (r), 
in the NL parametrization case and 

E l = E* + r„y 0 (r) ± irpMr) + eAo(r) 


dT v 

dpB 


PbVq + 


dT P p 3 
dp B 2 


bo + 


dT s 

dpB 


Ps3$0 


dT s 

a P s 
OPB 


(Jo 


in the DD case. The plus sign in the equations above 
has to be chosen for the proton and the minus sign for 
the neutron. Also, Vo, 4> 0 , &o and <5o are the time-like 
components of the meson fields. For the electron: 


E e = \Jk 2 + m 2 — e-4 0 (r); M* —> m e 
and for the neutrino: 

E v = k; M* -A 0. 

Finally, in order to take in to account the finite size of 
the nucleon, we have multiplied the constants gy ( gA ), 
or cy (ca), by the corresponding form factor. Using the 
parametrizations as explained in [2l|, we have taken a 
common structure factor given by: 


1 + 


4.97g 2 \ 

AM 2 J ’ 


(13) 


where q=\q\ = |pi -p 3 |- 
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